function [W, B, C, E, X, TOT,B2,Bex,Cex] = Decom(phi, sD, sD1,omega,omega1,sX1)
sDAGG = sum (omega .* sD);
sDAGG2 = sum (omega1 .* sD1);

sIAGG = 1 - sDAGG;
sIAGG2= 1 - sDAGG2;

I = ones(length(phi),1); % Identity matrix

si = I - sD; % import shares
si2 = I - sD1;

IM = (sD<1);
IM2 = (sD1<1);
X2 = (sX1>0);

CI = IM .* IM2;       % Continuing importers
NI = (I-IM) .* IM2;   % New Importers
OI = IM .* (I - IM2); % Old Importers

W = sum( CI .* omega .* (si2 - si)) / sIAGG * 100;
B = sum( CI .* (omega1-omega) .* si ) / sIAGG * 100;
C = sum( CI .* (omega1-omega) .* (si2-si) ) / sIAGG * 100;
E = sum( NI .* omega1 .* si2 ) / sIAGG * 100;  
X = sum( OI .* omega .* si ) / sIAGG* 100;
B2 = sum( CI .* (omega1-omega) .* si2 ) / sIAGG* 100;
TOT = (sIAGG2 - sIAGG) / sIAGG * 100;

Bex = sum( CI .* X2 .* (omega1-omega) .* si ) / sIAGG * 100;
Cex = sum( CI .* X2 .* (omega1-omega) .* (si2-si) ) / sIAGG * 100;

